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ABSTRACT 

Giant Molecular Clouds (GMCs) are massive, bounded, cool, dense regions containing 
mostly H 2 , but also Hi, CO, and other molecules. These clouds occupy < 1% of the 
galactic volume, but are a substantial part of the interstellar mass. They are irradiated 
by the high energy cosmic rays which are possibly modulated by the matter and magnetic 
fields within the clouds. The product of cosmic-ray flux and matter density is traced by 
the emission of high energy gamma-rays. In this paper, we consider a spherical cloud 
model and predict the gamma ray flux from several GMCs within 1 kpc of the sun which 
should be detectable by the EGRET instrument on GRO. 


I. INTRODUCTION 

Interaction of cosmic rays with matter in molecular clouds, containing mostly II 2 , but also 
III, CO, and other molecules, produces high energy gamma rays. These gamma rays pass nearly 
unattenuated through the galaxy and thus, are a good tracer of the product of the cosmic ray flux 
and the matter density at distant places in our galaxy. The density of Hi in these clouds can be 
observed directly. In cool, dense clouds, the Il 2 density can be inferred from the intensity of the 
CO, J = 1 — ♦ 0 transition. If the Hi and CO radio observations accurately represent the matter 
distribution in molecular clouds then the local cosmic ray flux at these clouds can be determined 
from the gamma ray flux. 

The gamma ray production rate as a function of the line of sight matter density, or source 
function, has been calculated for the Orion clouds by Bloemen, et al. (1984b). They used the CO 
data from Maddalena, et al. (1986), the Hi data derived from Heiles and Habing (1974, b < -10°) 
and Weaver and Williams (1973, b > -10°) and the COS-B gamma ray data, Mayer-Hasselwander, 
et.al. (1982), with a maximum-likelihood technique to derive the parameters in the relationship 


/ 7 = A-N{Ul) + B-W C o+C 


where A is the gamma ray production rate per H atom, B is the average A(H 2 )/VU C o ratio mul- 
tiplied by the gamma ray production rate per II 2 molecule and C is the isotropic galactic diffuse 
emission and instrumental background level. The CO and Hi data were smoothed to match the 
1 . 5 ° point spread function of COS-B. For 198° < / < 222° and -25° < 6 < -5° their results are: 


Energy Range 
A [ 7 sr -1 s -1 atom -1 ] 

B [7 cm -2 sr -1 K -1 (km/s) -1 ] 
C [7 cm -2 sr -1 s -1 ] 

A(H 2 )/VU C o = B/2A 
[mol cm -2 K -1 km -1 s] 


300 - 5000 [MeV] 
(0.52 ± 0.13) 10 -26 
(2.7 ± 1.0) 10 -6 
(2.0 ± 0.4) 10 -5 
( 2.6 ± 1 . 2 ) 10 20 


100 - 5000 [MeV] 
(1.70 ± 0.25) 10 -26 
( 10.1 ± 1 . 8 ) 10 -6 
(5.1 ±0.4) 10 -5 

(3.0 ± 0.7) 10 20 
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The value of /1(300 - 5000 MeV) found is consistent with the average emissivity value determined by 
Strong, et al. ( 1982) (0.59 ± ~ 10 %) • 10 -2G 7 s -1 sr -1 atom -1 using medium latitude galaxy counts 
and by Hlocmcn, et al. (198-la) (0.53 ± 0.14) ■ 10 _26 7 S “ 1 sr -1 atom -1 from the radial distribution 
of the gamma ray emissivity in the outer galaxy. 

The local cosmic ray source strength has most recently been derived by Strong, et al. (1988), 
by fitting the CO and Hi survey data, inverse compton emission, isotropic background and point 
source contribution to the COS-B diffuse galactic gamma ray emission. They did the fit for six 
gal actocentric rings, 2-4, 4-8, 8-10, 10 12, 12-15 and > 15 kpc (R @ = 10 kpc), and three energy 
ranges, 70-150, 150 300 and 300-5000 MeV. Their results, for 8 < R < 12 kpc, are 

{ 1 .02 ± 0.10 X 10 -2G 7 sr -1 s -1 atom 1 (70 < E < 150 MeV) 

0.65 ± 0.06 x 10 -26 (150 < E < 300 MeV) 

0.62 ± 0.06 x 10 - 26 ( 300 < E < 5000 MeV) 

Most source strengths arc quoted for the energy range E> 100 MeV. Assuming a power law spectra 
with index — 2 . 1 , the above results can be expressed in the form 


% 


3.7 x 10 5 7 mg 1 s 1 sr 1 
1.0 x 10 -5 


(20 < E < 100 MeV) 
(E > 100 MeV) 


We will use the this source function and a spherically symmetric CMC model with radially 
dependent density to predict the gamma ray emission of GMCs which will be observable by EGRET. 


II. SPHERICAL CLOUD MODEL 

Following Solomon et al. (1987), to model the gamma ray emission from GMC’s, we consider 
the simplified case of a spherically symmetric cloud with radial density of the form 


/>(r)[g c m 3 ] = po ■ (IWr)° 


where Rq is the cloud radius. The line of sight column density is given by 


^( r )[g 



p(r = \/r 2 + x l ) di 


where r is the radial distance between the cloud center and the line of sight, or “impact parameter”. 
Figure 1 shows the line of sight column density for a = 0, 1 and 2 normalized for clouds with 
M = 10 6 Mq and R = 50 pc as a function of r. A similar model for cloud density was used by 
MacLaren, et al. (1988) who normalized the density at r/Ro = 0.55. 

Analysis of the CO data by MacLaren, et al. (1988), shown in figure 2 , suggests that a is close 
to 1 , but it may be almost 2 for some clouds. This can also be seen by comparing the a = 1 column 
density with the W c o contours of the Orion A and B cloud peaks in figures 6 and 7 of Maddalena, 
et al. (1986). The peak at (a = 5 h 44 m , 6 = 0°), for example, has a fairly linear density fall-off with 
radial distance from the peak and a long tail. We include a = 0 and 2 as limiting cases. 


III. PREDICTED EMISSION OBSERVED BY EGRET 

The diffuse galactic gamma ray spectrum above about 10 MeV is composed of an E -2,5 Brehm- 
strahlung component and the 7 r° decay spectrum centered at about 70 MeV. We will consider here 
gamma rays with energies above about 100 MeV, where the ir° spectrum begins to dominate. The 
effective area of EGRET in this energy region is about 1460cm 2 . 
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Figure 1. Column density as a function of impact parameter for a cloud with 
density given by p — po -(Ro/r) 01 , for ot = 0, 1 and 2, normalized to the same 
total mass, M = 10 6 M 0 , and radius, R = 50 pc. 
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Figure 2. Column density of hydrogen (arbitrary units) derived from CO 
data vs. normalized cloud radius for several local molecular clouds. The 
solid curves are column densities for spherical clouds with density given by 
p = p 0 (Ro/r) a , for a = 0, 1 and 2, normalized at r/R 0 = 0.55. Orion(A+B), 
(+); Taurus, (o); Perseus, (x and □). From MacLaren et al. (1988). 


45 





ORIGINAL PAGE IS 
OF POOR QUALITY 



Figure 3. Number of gamma rays observed from a hypothetical GMC with 
cloud density given by p = p 0 * (Ro/r)°y for a = 0,1, and 2, M = 10 6 M 0 , 

Re — 20 pc at D = 500 pc (left) and 1000 pc (right). The total counts for 
T 0 bs = 6 X 10 5 s are given per 0.5° ring. 

To match the EGRET instrument angular resolution of about 0.5° for gamma ray energies 
above 100 MeV, it is appropriate to sum the total number of gamma-rays observed from a cloud 
into rings of width 0.5°. Thus, the total line of sight mass of a ring with inner radius r x and outer 
radius 7*2 is given by 

S( r i , r 2 )[g] = 2 tt f r • cr(r) dr 
Jr x 

If all the sources of gamma ray production in the cloud can be expressed in terms of a single source 
function, such as that derived by Strong, et al. (1988), then the approximate gamma ray flux from 
the ring is 

Counts/ring[7 sr -1 s -1 ] = / • q 0 • £(77, r 2 ) 

where q 0 [7 g _1 sr -1 s -1 ] is the local cosmic ray source function and / is the relative cosmic ray 
intensity at the cloud. The calculation of the number of observed gamma ray from a GMC requires 
knowledge of the total mass of the cloud and the distance to the cloud in addition to the sensitivity 
or effective area of the detector. 


CountSobs/ring = A e n • T 0 b s • 4n/4irD 2 - Counts/ring 

where A e ff [cm 2 ] is the energy dependent effective area, T 0 b s [s] is the observation time and D [cm] 
is the distance to the cloud. 

In figure 3 the number of counts, E > 100 MeV, per 0.5° ring are shown for a hypothetical 
10 6 M© cloud with 50 pc radius immersed in cosmic rays of the local density (/ = 1) for an 
observation time of 6 X 10 5 s (effective time for a two week observation) at a distance of 500 and 
1000 pc respectively. The diffuse galactic and extra-galactic gamma ray emission has not been 
included. Significant differences in the count profile can be seen for the different values of a. For 
a - 0 the cloud appears very much as a ring with a sharp edge. For a = 1 the ring structure is 
still visible, but the edge is softer. For a - 2 the cloud appears as a bright spot with a soft edge. 
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We have calculated the number of gamma rays which EGRET can expect to observe in a 
typical two week observation for the clouds listed in table 1. These clouds, selected from the work 
by Dame et al. (1987), are more than 25 pc above or below the galactic plane, within 1 kpc of the 
sun and have masses greater than 1 x 10 4 M 0 , see figure 4 . 


1 ° ■ 
mtn 


Cygnus OB 7 

87.0 

Cepheus 

100.0 

Taurus 

163.0 

Orion B 

202.5 

Orion A 

208.5 

Chamaeleon 

295.0 

From Dame et al. 

(1987) 


TABLE 1 


1° 

‘max 

b°min 

b° 

u max 

99.0 

-3.0 

8.0 

120.0 

11.0 

22.0 

178.0 

-22.0 

-9.5 

208.0 

-21.0 

-6.0 

218.0 

-21.0 

-14.5 

305.0 

-20.0 

-12.0 


D[pc] M[1O 5 M 0 ] 


SOO 

7.5 

450 

1.9 

140 

0.3 

500 

1.7 

500 

1.6 

215 

0.1 



Figure 4. The distribution of GMCs with mass greater than 1 X 10 4 M 0 
within 1 kpc of the sun. Adapted from Dame et al. (1987). 
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Combining the emission data of Strong et al. (1988) with the mass, radius and distance data 
in Table 1, we can calculate the high energy (E > 100 MeV) gamma ray emission from these clouds. 
Figure 5 a b and c, shows the estimated number of gamma rays (E > 100 MeV which will be 
observed from each of these clouds during a two week observation (6 X 10 5 s) summed into 0.5° 
circular bins. In each figure, the diagonal dotted line is the estimated number of counts from 
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the diffuse emission for the same observation. The diffuse gamma ray emission (E > 100 MeV 
was estimated using the data given by Fichtel et al. (1978). The appropriate longitude data was 
averaged over the width of the cloud in latitude. 

We now consider the statistical significance of these observations in two ways. The significance 
of the gamma ray observed in a ring over the background reaches a broad maximum at about 30% 
of the cloud radius and falls to zero at the cloud edge, as defined by the model. The significance of 
the total number of gamma* rays observed from a cloud over the background reaches a. maximum 
at about 82% of the cloud radius and then starts to fall as the background begins to dominate the 
signal. These values arc model dependent, however, it is clear that knowledge of the cloud radii will 
be important in the analysis of a wide field of view gamma ray observation to determine the gamma 
rays which are produced in molecular clouds. Table 2 gives the modeled detection significance (in 
standard deviations) of the ring of maximum significance and of the total cloud for the six clouds 
listed in table 1. 


TABLE 2 

Significance of 
Maximum Ring Total Cloud 


Cygnus OB7 

2.2 

5.5 

Cep heus 

14.0 

39.G 

Taurus 

23.5 

64.0 

Orion B 

21.4 

47.1 

Orion A 

24.4 

50.7 

Chamacleon 

6.9 

15.1 


(in standard deviations) 


IV. CONCLUSION 

Clouds with masses greater than about 10 5 M 0 and closer than about 1 kpc should be detectible 
with the EGRET instrument in a two week observation. In addition, local and cloud-to-cloud 
variations in the product of the column density and the cosmic ray density should be resolvable 
for many of these clouds. We expect to also be able to detect some clouds to distances of a several 
kpc, and should be able to use these clouds as tracers of the cosmic ray density in more distant 
parts of the galaxy. 
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Figure 5. a Estimated gamma ray counts (E > lOOMeV) from GMCs in Ce- 
phus (left) and Taurus (right) summed into 0.5° circles for a two week ob- 
servation. The dotted line is the estimated diffuse emission for the same 
observations. 



Figure 5.b Estimated gamma ray counts (E > lOOMeV) from GMCs in Cha- 
maeleon (left) and Cygnus (right) summed into 0.5° circles for a two week 
observation. The dotted line is the estimated diffuse emission for the same 


observations. 
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Figure 5.c Estimated gamma ray counts (E > lOOMeV) from OrionJV (left) 
and Orion.B (right) summed into 0.5° circles for a two week observation. 
The dotted line is the estimated diffuse emission for the same observations. 
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